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The PHENIX experiment has measured electrons and positrons at midrapidity from the decays 
of hadrons containing charm and bottom quarks produced in d+Au and p+p collisions at y/s NN = 
200 GeV in the transverse-momentum range 0.85<py<8.5 GeV/c. In central d+Au collisions, the 
nuclear modification factor RdA at 1.5 < pr < 5 GeV/c displays evidence of enhancement of these 
electrons, relative to those produced in p+p collisions, and shows that the mass-dependent Cronin 
enhancement observed at the Relativistic Heavy Ion Collider extends to the heavy D meson family. 
A comparison with the neutral-pion data suggests that the difference in cold-nuclear-matter effects 
on light- and heavy-flavor mesons could contribute to the observed differences between the 7T° and 
heavy- flavor-electron nuclear modification factor -Raa- 

PACS numbers: 25.75.Cj 



The experimental collaborations at the Relativistic 
Heavy Ion Collider (RHIC) have established that a hot, 
dense medium with partonic degrees of freedom is formed 
in Au+Au collisions at y/s^ = 200 GeV The tem- 

perature achieved in this medium, as inferred from direct- 
photon measurements, is well over the threshold expected 
from latticc-quantum-chromodynamics (QCD) calcula- 
tions to enable deconfinement and create the quark gluon 
plasma [f|. Studies of the interactions of heavy quarks 
with this matter are of particular interest. Since charm 
and bottom quarks are dominantly produced by gluon 
fusion in the early stages of the collision, they experience 
the complete evolution of the system. The heavy-quark- 
production baseline in p+p collisions at y/s NN = 200 GeV 
is consistent with fixed order plus next-to-leading-log per- 
turbativc QCD calculations within uncertainties [6j- In 
central Au+Au collisions, suppression of electrons from 
the decays of hadrons containing heavy quarks has been 
measured relative to the yield in p+p scaled by the num- 
ber of nucleon-nucleon collisions, iV co ii, suggesting that 
heavy quarks lose a significant amount of their initial 
energy [7[. The positive elliptic flow amplitude of these 
decay electrons implies that heavy quarks flow along with 
the light partons that compose the bulk of the medium. 
When considered together, the suppression and elliptic 
flow of these quarks are in qualitative agreement with 
calculations based on Langevin transport models that 
imply a viscosity to entropy density ratio close to the 
conjectured quantum lower bound of l/4ir (3-HoT|. 

A full understanding of these phenomena requires mea- 
surements of cold nuclear matter (CNM) effects, which 
are believed to be present in Au+Au collisions but are dif- 
ficult to distinguish experimentally from effects due to in- 
teractions with the hot medium. Initial state effects, such 
as modification of the parton distribution functions in the 
nucleus, momentum broadening due to parton scattering 
in the nucleus, and CNM energy loss will affect heavy- 
quark production rates and cannot be accounted for with 
a reference from p+p data 11 - 13] ■ It is therefore neces- 
sary to study p+Au (or d+Au) collisions, where a hot 



nuclear medium is not expected to form, to isolate these 
nuclear effects. 

To this end, a vigorous experimental effort to quantify 
CNM effects is underway at RHIC. A mass-dependent 
Cronin enhancement has been observed for tt, K, and p 
production , where the pt spectra of these hadrons 

in d+Au collisions are hardened with respect to p+p. 
While overall J/ip production is suppressed in d+Au col- 
lisions, a broadening of the pt spectrum is also observed 



16l Il7| . The relative strengths and centrality depen- 



dence of initial-state effects and breakup in the cold nu- 
clear medium that contribute to these phenomena are 
not known. The study of mesons containing open heavy 
flavor can help disentangle these coexisting effects. This 
Letter presents measurements of pr spectra and the nu- 
clear modification factor (RdA) of electrons and positrons 
from the decays of hadrons containing charm and bot- 
tom quarks (e^ F ) produced in d+Au collisions at ^/s NN 
— 200 GeV. When combined with heavy-quark measure- 
ments from p+p and Au+Au, this analysis provides a 
detailed study of the production of heavy quarks, the ef- 
fects of production in a nucleus, and the dynamics of the 
hot nuclear medium. 

The PHENIX experiment ^ sampled 80 nb" 1 of in- 
tegrated luminosity during the 2008 d+Au run at RHIC, 
a factor of 30 increase over the 2003 d+Au data set. 
The minimum bias (MB) trigger and event centrality 
are obtained from two beam-beam counters located at 
3.1 < | r) | < 3.9 in pseudorapidity. The charge generated 
in the beam-beam counter facing the incoming Au nu- 
cleus is divided into four categories covering the 0-20%, 
20-40%, 40-60%, and 60-88% most central collisions. As 
the MB-trigger efficiency is 88+4% of the total d+Au in- 
elastic cross section, a correction factor is applied to the 
yield measured in the MB-triggered data sample to give 
a nonbiased sample, covering 100% of the d+Au collision 
centrality. 

This analysis considers electrons and positrons identi- 
fied in the two PHENIX central arm spectrometers. Each 
arm covers an azimuthal angle A</> = n/2 and a pseudo- 
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rapidity range \r]\ < 0.35, and uses layers of multiwire 
proportional chambers and pad chambers for charged 
particle tracking. Ring-imaging Cerenkov (RICH) coun- 
ters and electromagnetic calorimeters (EMCal) provide 
electron-identification and hadron-rejection capabilities. 
A coincidence of the MB trigger and a RICH hit matched 
with an energy deposit of at least 600 or 800 MeV in 
the EMCal functions as an electron trigger. At px = 5 
GeV/c, charged pions begin to radiate in the RICH, but 
matching requirements between the track's energy de- 
posit in the EMCal and reconstructed momentum ef- 
fectively eliminate hadron contamination out to px = 8 
GeV/c. Above this, hadronic contamination accounts for 
20 ± 10% of the signal, and is subtracted. A full GEANT 
simulation of the PHENIX detector is used to correct 
for the incomplete azimuthal acceptance and electron- 
identification efficiency of the central-arm detectors. 

Most of the electrons produced in collisions at RHIC 
come not from heavy-flavor decays, but from the neutral- 
pion Dalitz decay, ir° — » ~fe + e~. The r\ Dalitz decay 
contributes about 10% of the electron background for 
1 < pt < 9 GeV/c. Other hadron decays (77', p, uj, cj>, T) 
add to the background at the few percent level. In- 
ternal and external conversions of direct photons, while 
negligible at px < 2 GeV/c, are significant sources of 
electrons at high momentum. Electrons from the decay 
J/ip — » e + e~ are a significant source of background at in- 
termediate pt, and constitute a maximum of about 25% 
of the total electron background at px — 5 GeV/c. Con- 
versions of photons from hadron decays are significant 
at all momenta, however, the low material design of the 
PHENIX detector ensures that the number of these con- 
version electrons is less than half of that from neutral- 
pion Dalitz decay. In addition, electrons produced at 
displaced vertices from the K e 3 decays of K mesons are 
misreconstructed by the PHENIX tracking algorithm and 
contribute about 3% of the total background at pt = 0.85 
GeV/c, but quickly fall off to less than 1% at px = 1-5 
GeV/c. 

Two independent methods are used to isolate the con- 
tribution of heavy flavor electrons. The cocktail method 
uses a Monte Carlo hadron decay generator to calcu- 
late the electron background from each relevant hadron 
species. The parametrization of the neutral-pion px 
spectrum is determined by a modified Hagedorn fit to 
pion data obtained from earlier measurements in d+Au 
[3, Ell ■ The shape of the px spectra of the other mesons 
is determined by wit scaling the pion fit, that is, the vari- 
able substitution px — > mx — \Jpt 2 + (^mcson ~ m ^o), 
and their normalization is set to world av erag es of the 
ratio of mcson/7r° at high momentum [l9|, |20| Direct- 
photon contributions are estimated by scaling the mea- 
sured direct- photon yield in p+p by A co n [2l| • The num- 
ber of conversion electrons is found by a full GEANT sim- 
ulation of the PHENIX detector material, and a simi- 



lar simulation, in conjunction with the actual PHENIX 
tracking algorithm, is used to estimate the K e $ decay 
background. Contributions from J/ip decays are found 
by parameterizing the measured J / ip spectrum from [l6| 
for each centrality, for d+Au, and from [22[ for p+p. The 
small background due to T decays and the Drell-Yan 
process are taken from [23[, and scaled by N co u for each 
centrality. The sum of these background sources is then 
subtracted from the inclusive electron measurement to 
give the heavy flavor contribution. 

The second method of signal extraction is based on the 
fact that the vast majority of the background electrons 
are "photonic" in nature, i.e. they originate from either a 
real photon (the conversion electrons) or a virtual photon 
(the electrons from Dalitz decays) , while signal electrons 
are nonphotonic. The inclusive yield of electrons in the 
standard detector configuration can be parametrized as 

^standard _ ^7 _|_ ^ynorry ^ 

where TV 7 (TV 11 ™ 7 ) represents the photonic (nonphotonic) 
electron yield. The addition of extra material (the "con- 
verter" , a sheet of brass 1.68% of a radiation length thick, 
wrapped around the beam pipe) into the PHENIX aper- 
ture increases the photonic component by a factor i? 7 , 
but attenuates the signal by an amount (1-e), giving a 
total yield 

^converter = + (1 - e)7V n ° n7 (2) 

By modeling the converter material in simulation, the 
factors R 7 and e are determined to be 2.32 ± 2.7% (with 
a slight px dependence), and 0.021 ± 25%, respectively. 
The inclusive yields 7Vf andard and N™ nveItel are mea- 
sured by the PHENIX spectrometer, so a simultaneous 
solution of Eqs. (1) and (2) gives the quantity of inter- 
est iV non7 . The nonphotonic background sources, namely 
K e 3 decays and the dielectron decays of the p, ui, <p, J/tp, 
and T contribute about 10% of the total background at 
px < 1 GeV/c, and are subtracted following the cocktail 
method described above. The converter method provides 
a robust but statistics-limited determination of the pho- 
tonic background. Since the converter material creates 
an undesirable background for other measurements, only 
3% of the tf+Au data recorded by PHENIX in 2008 was 
taken with the converter installed. 

A crucial cross-check of this measurement's accuracy 
is the consistency of these two independent background 
determination methods. A comparison of the photonic 
components of the cocktail (Dalitz decay electrons, con- 
versions, and direct photons) to the photonic-electron 
signal extracted by the converter method shows agree- 
ment within 8% for all centralities (see inset of Fig. [TJ. 
Since the converter method gives a direct measurement of 
the photonic background, while the cocktail is a calcula- 
tion that relies on simulation, the photonic components 
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FIG. 2: (color online) The nuclear modification factor, RdA, 
for electrons from open heavy flavor decays, for the (a) most 
central and (b) most peripheral centrality bins. 



FIG. 1: (color online) Electrons from heavy flavor decays, 
separated by centrality. The lines represent a fit to the previ- 
ous p+p result [53], scaled by AT co ii. The inset shows the ratio 
of photonic background electrons determined by the converter 
and cocktail methods for Minimum Bias d+Au collisions, with 
error bars (boxes) that represent the statistical uncertainty on 
the converter data (systematic uncertainty on the photonic- 
electron cocktail). See text for details on uncertainties. 



of the cocktail are scaled to match the converter data 
in each centrality by factors ranging from 0.92 to 1.01. 
Detailed descriptions of these methods can be found in 
& 

Figure [T] shows the px spectrum of electrons from 
open heavy flavor decays for each d+Au centrality bin, 
and for p+p collisions that were measured during the 
same RHIC Run period with identical techniques. The 
heavy flavor electron yield is determined by the cocktail 
method, with photonic components scaled to match the 
converter data. The statistical (systematic) uncertain- 
ties are shown as bars (boxes) around the central values. 
The boxes contain the uncertainties in the solid angle 
correction, electron-identification efficiency, and trigger- 
bias correction. Added in quadrature with those is the 
uncertainty from the cocktail subtraction. The lines are 
a FONLL spectral shape 24| fitted to a previous p+p 
heavy- flavor electron measurement [23|, scaled by N co u 



for each centrality. The p+p data presented here are in 
good agreement with our previous p+p results, however, 
the statistical uncertainties on the new data are ~ 2x 
larger. Fitting a constant to the ratio of the new data 
to the old yields a value of 0.97 ± 0.02, with j^/cudi = 
20.3/26. The fact that the 2008 p+p data agree with the 
previous p+p data provides an important cross-check on 
the methods used to extract the 2008 d+Au e^ F spectra. 

Due to changes in the detector configuration that re- 
sulted in increased photon conversion background at low 
Pt, the signal to background at low px is not as good as 
it was in previous measurements. Coupled with the fact 
that ~ 90 % of the electrons from charmed hadron de- 
cays fall below pt — 0.8 GeV/c, where the present data 
cut off, this means that the data do not place meaningful 
constraints on the total charm production cross section. 

The d+Au electron spectra are directly compared to 
the p+p reference data by computing 



RdA 



dNlJdp T 



(N con ) x dN'/dpi 



(3) 



for each centrality. Figure [2] shows RdA as a function 
of pt for the most-peripheral and most-central central- 
ity bins. As in Fig. [IJ the statistical (systematic) un- 
certainties are represented by bars (boxes). For points at 
Pt < 1.6 GeV/c, RdA is found by dividing point- by-point 
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FIG. 3: (color online) The nuclear modification factors RdA 
and Raa for minimum bias d+Au and Au+Au collisions, for 
the 7t° and e HF . The two boxes on the right side of the 
plot represent the global uncertainties in the d+Au (left) and 
Au+Au (right) values of iV C oii- An additional common global 
scaling uncertainty of 9.7% on RdA and Raa from the p+p 
reference data is omitted for clarity. 

the d+Au yield by the p+p yield from [23]. At higher 
transverse momentum, where the p+p heavy-flavor elec- 
tron spectrum is consistent with a shape from pQCD, 
a fit to the spectral shape from the [24| calculations is 
used to represent the p+p yield. The statistical uncer- 
tainty on the fit is included as a systematic uncertainty 
on the shape of RdA by adding it in quadrature with the 
systematic uncertainties on the electron background sub- 
traction and solid angle and efficiency corrections. The 
global scaling uncertainty from the uncertainty in N co ii 
and the total sampled p+p luminosity is given by a box 
on the right. Note that the 2008 p+p data shown in Fig.Q] 
could be used for the denominator of RdA, however, the 
use of the more precise data from [23| gives a smaller 
uncertainty on RdA- 

The central RdA shows an enhancement out to pt ~ 5 
GeV/c, and implies that the suppression of heavy flavor 
electrons in central Au+Au collisions at RHIC is not an 
initial state CNM effect, but rather is due to the hot nu- 
clear medium. The peripheral nuclear modification fac- 
tor also shows some evidence of an enhancement, which 
is to be expected since even the most peripheral central- 
ity bin in d+Au samples a significant nuclear thickness. 
Although the techniques used here do not allow sepa- 
ration of electrons from charm and bottom decays from 
each other, measurements from p+p show that pr = 5 
GeV/c is near the transition point where contributions 
from bottom quarks begin to dominate over charm (jjj . 
Since the total charm cross section is expected to scale 
with iV co ri, this enhancement below 5 GeV/c suggests a 
Pt broadening of the D spectral shape, with a mass de- 
pendence that roughly follows the previously observed 
trend in the 7r, K, and p families. The B spectrum may 



also be modified, however, the uncertainties on the data 
and on the relative D and B contributions to the electron 
spectra preclude a precise determination of any effects. 

The effects of cold nuclear matter are expected to be 
present in the initial state of A+A collisions, however, 
this CNM enhancement is convolved with the suppress- 
ing effects of hot nuclear matter. Figure [3] shows RdA and 
Raa for e^ F and the neutral pion, for which only small 
CNM effects are observed [II El- Above p T ^ 5 GeV/c, 
where the CNM effects on both species are small, their 
Raa values are consistent within uncertainties. How- 
ever, in the range where CNM enhancement is large for 
e^ F and small on 7r°, the corresponding e^ F i?AA values 
are consistently above the n° values. This could suggest 
that the difference in the initial state cold nuclear matter 
effects due to the mass-dependent Cronin enhancement 
is reflected in the final state spectra of these particles 
in Au+Au collisions, although alternate explanations in- 
volving mass-dependent partonic energy loss in the hot 
medium are not ruled out. 

In summary, we have observed an enhancement of elec- 
trons from heavy-flavor decays produced in central d+Au 
collisions at y/s NN = 200 GeV. The previously observed 
suppression of these electrons in central Au+Au collisions 
is therefore attributed to hot-nuclear-matter effects. We 
find that the ir° and ej| F nuclear modification factors 
Raa are consistent within uncertainties in the px range 
where CNM effects on both species are small. In the 
range where CNM enhancement of e^ F is significant in 
d+Au, these effects may also be apparent in the Au+Au 
data. 
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